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■ Abstract 

, We present a model for the structure of the pion. Based on ideas of a 

^ I recently developed statistical model of the nucleon, we assume the pion to be 

a gas of partons. The finite-size corrections (FSC) are incorporated through 
, two parameters. Using the same two FSC parameters for the proton and 

I pion we reproduce quantitatively the data on ir~N fj,'^fi~X Drell-Yan 

'sj" ■ production and valence quark distribution of the pion. 
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1 Introduction 



Recently, a statistical model for the structure of the nucleon was proposed The 
nucleon was described as a non-interacting gas of valence quarks, sea quarks and glu- 
ons in equilibrium. This picture was then improved by adding finite-size corrections 
(FSC) to the expression for the parton density function (PDF). These finite size cor- 
rections take into account the fact that the gas is enclosed in a finite volume. The 
model was successful in reproducing a large body of polarized and unpolarized nu- 
cleon structure function data. Since the ideas from statistical mechanics worked well 
inside one hadron, namely the nucleon, it is quite motivating to test these ideas in 
other hadrons as well. However, compared to the vast amount of precise data avail- 
able for the nucleon, much less is available for others. Among hadrons other than the 
nucleon, the pion is the most widely explored. There exist data on vr^A^ ^^fi^X 
Drell-Yan production |0, |^ and prompt photon production in n^p 7X which 
have been used by some groups f^, § to obtain parameterizations for PDFs of the 
pion. The Drell-Yan data constrains the shape of the pion valence densities and the 
prompt photon data constrains the pionic gluon distribution in the large-x region. 
These data are, however, not sufficient to fix the gluon and sea in the pion uniquely. 
The number of parameters appearing in the expressions for the PDFs obtained from 
such global fits are very large and with no physical meaning. Hence it is worth 
investigating the scope of a statistical model which has fewer parameters based on 
a physical interpretation. 

In this Letter, we present a statistical model for the parton distributions in the 
pion. The expressions for the pion PDFs are written assuming the pion to be a gas 
of massless partons. These are then improved to include FSC as in the case of the 
nucleon. In principle the two parameters associated with the FSC can be determined 
theoretically [0 . However, these parameters are sensitive to the equation of motion 
employed, the boundary conditions on the wave function, shape of the enclosure and 
details such as whether the particles are strictly massless. In ref . , in the absence 
of a complete understanding of the nucleon as a QCD bound state, a practical 
approach was taken and the two parameters for the nucleon were determined by 
fitting unpolarized structure function F2{x, Q'^) data at one value of Q"^. In the case 
of the pion, due to scarcity of precise experimental data giving direct information 
on the pion PDFs, we prefer to use as a first guess, the same parameters for the pion 
as in the case of the proton [Q. Interestingly, with these parameters we get good 
quantitative agreement with the available data on pion valence quark distribution 
and 7r~iV — Drell-Yan production. The agreement with data is nearly as 
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good as that obtained by the existing parameterizations 0. 



2 The model 

We visuahse the pion to be made up of a gas of massless partons in equhbrium at 
temperature T in a spherical volume V with radius R. The parton number density 
dn'^/dx in the infinite momentum frame (IMF) and the density dn/dE in the pion 
rest frame are related to each other by, 

dv} M'^x /-A^/a dE dn 

(1) 



dx 2 JxM/2 E"^ dE 

where the superscript i refers to the IMF, M is the pion mass and E is the parton 
energy in the pion rest frame fl], Using the standard procedure in statistical 
mechanics to introduce the effects of the finite size of an enclosure, in the expression 
for the density of states P| we write dn/dE for the pion as, 

dn/dE = g f{E) iyE'^/l-n^ + aWE + (2) 

where g is the spin-color degeneracy factor, jiE) is the usual Fermi or Bose distri- 
bution function f{E) = {exp[{E — fi)/T] ± 1}^^, V is the pion volume and R is the 



radius of a sphere with volume V. We take R = yS/Sp = 1.07 fm, where p is the 
root-mean-square radius of the pion |[T^. The three terms in (2) are the volume, 



surface and curvature terms, respectively. The coefficients a and b are, in principle, 
the free parameters of the model. For reasons mentioned in the introduction, though 
a and b can be determined theoretically, we prefer to choose these values to be the 
same as those determined phenomenologically for the proton in ref. []T| . Theoretically 
this makes sense, as the values a = —0.376 and b = 0.504 obtained from the fitting 
procedure in ref. are close to the values (a = —1/2 and b = 3/27r) determined 
theoretically by Morse and Ingard 0. We will come back to this point at the end 
of the paper. 

To determine the values of the temperature and chemical potentials appearing 
in the distribution function f{E), we consider the fact that any model of the PDFs 
for the pion has to obey the number and momentum constraints, the former being 
motivated by the constituent quark picture for mesons at a low scale (this scale, 
being our input scale, should not be to too small as compared with Aqcd to allow 
for the evolution of the partons densities, hence we choose it to be = where 
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mjv is the mass of the nucleoli). If na(Q-) denotes the number of quarks (antiquaries) 
of flavour a, then in the case of a tt^ for example, the constraints can be given as, 

Uu - Uu ^ 1 (3) 

- = 1 (4) 

Us - ris = (5) 

{momentum fractions) — 1. (6) 

The numbers na(a) in eqs.(3-5) are obtained from eqs.(l) and (2) by integrating the 
appropriate dn'^/dx over x. The momentum fractions are obtained by integrating 
xdn^ /dx over x. The temperature T and chemical potentials (/x) appearing in eqs.(3- 
6) are then not free parameters. For a given set of parameters a and 6, they are 
determined uniquely by solving the four coupled non- linear eqs.(3-6). 

Using the values of a and h mentioned above, we obtain T = 22 MeV, = 
—/lu = 162 MeV and /x^ = in the case of a negatively charged pion. The left and 
right hand sides of eqs.(3-6) with these values of T and /x, agree with each other 
up to an accuracy of one part in 10^. The chemical potentials for antiquarks are 
determined by using the relation /jq = —jJLq. The model, described above, fixes the 
parton distributions in a pion uniquely at an input scale Qq. It remains to put it 
to test and compare the predictions with available experimental data. In order to 
compare with data, the parton distributions are evolved to different Q"^ values using 
the standard DGLAP evolution equations at next-to-leading order, starting from 
the input scale Qq = m^. 



3 Results and discussion 

To start with, we compare the vr^ valence quark distribution calculated within the 
statistical model with the available data at = 30 GeV^. In Fig. la we show the 
valence structure function extracted from n'W n^n~X Drell-Yan data by the 
E615 collaboration. The sohd curve corresponding to our calculation of xv^^ — x{u — 
u) done within the statistical model, shows reasonably good agreement with the E615 
data. Our values of xVj^ are, however, higher than those due to the SMRS (dashed 
line) and GRS (dash-dotted line) parameterizations. Both these parameterizations 
have determined their valence densities by making a global fit to the irN Drell-Yan 
data. To further compare our valence densities with other calculations, we calculate 
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the first two moments of tlie pion valence quark distributions given by, 



< xVt, > = / xvT,{x)dx (7) 
Jo 

< x^f^r >= / x"^ Vt,{x) dx. (8) 
Jo 

In Fig. lb we plot tlie dependence of the pion valence moments. The solid 
curve corresponds to our calculation and the SMRS result is shown by the hashed 
region. The points with error bars are the results from a lattice QCD calculation 
[ig performed at = 49 GeV^. 

Next, we put the statistical model to test with the 7r~N — > ^~^fi~X Drell-Yan 
data. We calculate the double-differential cross-section d'^a/dxpd^/r where xp = 
and ^/T = XttXn = M'^/s. At leading order, x^r and xn are the Bjorken 
X variables of the pion and target nucleon respectively. M is the invariant mass 
of the muon pair and ^/s is the centre of mass energy. We calculate the Drell- 
Yan cross-sections using the full next-to-leading order expressions § and with the 
choice = M^. The proton PDFs required for the calculation of the 7r~N Drell- 
Yan cross-sections are also calculated within the statistical model. Both the pion 
and nucleon PDFs are calculated using the same values of the FSC parameters a 
and b. A detailed description of the evaluation of the nucleon PDFs can be found 
in ref. . To take into account uncertainties due to normalization and higher order 
QCD corrections, we multiply our cross-sections by the standard i^'-factor. We 
analyse two different sets of data obtained from a vr" beam incident on a tungsten 
target, tt^W — > fi~^fi~X, by the NAlO and E615 collaboration. To make a correct 
comparison of our calculated cross-sections for n^N ^ fi^fi'X with the 7i~W data, 
we multiply our cross-sections with a factor R = —0.55a; at + 1.1. This form of the 



correction factor R is consistent with the observed values [|I^ of R as shown in 
ref.0. In Fig.2 we compare our calculated results with the NAlO data taken at tt~ 
beam momenta of 194 and 286 GeV/c. With K' factors of 1.02 and 1.06 at beam 
momenta 194 and 286 GeV/c respectively, we find reasonably good agreement with 
data at different values of y/r. Our K' values are close to those in ref.p which 
range between 1.1 and 1.4 depending on the type of fit. A discussion on the K' 
factor can be found in ref.[0] where it was also noted that it should be fairly close 



to unity. In Fig.2 we also show the dependence of the cross sections on the choice of 
the parameters a and b. The solid curves correspond to a = —0.376 and b = 0.504 
which are the values used throughout this work. Varying the parameters arbitrarily 
we found that the agreement with data reduces as we move away from the values of 
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a and b quoted above. To demonstrate this fact, we plot the cross sections at beam 
momentum 286 GeV/c for two different sets of a and b in Fig. 2. The dashed curves 
correspond to the calculations with a = —0.2, b = 0.3 and the dot-dashed curves 
correspond to those with a = —0.6, b = 0.7. 

In Fig. 3 we show our calculations in comparison with the E615 data at Pl = 252 
GeV/c. The K' factor in this case is 0.85. It differs from the K' obtained by 
us for the NAlO data due to normalization errors in the data. Since the ratio 
K' {N AlO) / K' {E615) in this work is very similar to that obtained in refs.0, 0, we 
do not worry about the fact that K'{E615) < 1. As mentioned above it could be 
due to errors in normalization. The agreement of our results with data is again 
quite good. We are, however, not able to reproduce data at very high or low values 
of y/r very well. Although the reason for this disagreement is not clear to us, it 
certainly is not a short-coming of the statistical model. We say so because this kind 
of disagreement was also found in ref. |Q. They included the data only in a certain 
range of i/r for fitting as they were unable to reproduce the data at very high and 
low values of a/t. 

Finally in Fig. 4 we compare the valence, sea and gluon distributions of the 
pion within the statistical model with those obtained by the SMRS and GRS pa- 
rameterizations. Since the pion sea distributions are not well determined by data, 
the SMRS parameterization made different fits by varying the fraction of the pion 
momentum carried by the sea between 10% and 20%. Interestingly, we find that the 
sea calculated within the statistical model carries a fraction roughly equal to 15% 
of the pion momentum which is like the SMRS best fit value. 

To summarize, we can say that we have presented a simple physical model for the 
pion, using ideas from statistical mechanics. The only parameters in the description 
of the pionic parton densities are the two 'finite-size correction' parameters. Con- 
sidering the fact that the only two parameters appearing in the model for the pion 
have not been determined by any fit to the pionic data, the success of the model in 
reproducing data on the pion is surprising. Since we use the same parameters for 
the proton and pion (for reasons explained in the Introduction), we have in a sense, 
a unified model for the proton and pion. We tried different sets of the values a and 
b for the pionic partons and having the same values as that for the proton is an 
outcome of our work rather than an assumption. At this stage, however, it would 
be too early to speculate about the significance of this result. It could be that given 
a wealth of precise enough data on pionic parton distributions, the 'true' values of 
a and b determined by a fit, will differ from the ones presently used. However, it 
is clear from our results that the difference will not be drastic. It seems therefore 
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that the values for a and h are always close to the theoretical values, a = —1/2 and 
b = 3/2-7? determined in ref. 0. This could partly be the reason why the statis- 
tical model works so well for the proton as well as for the pion using the same a 
and b. Further efforts to understand why the ideas from statistical mechanics work 
inside hadrons could be worthwhile and could improve our understanding of hadron 
structure. In future we intend to understand this point in more detail, by using 
for comparison a larger set of data (i.e. prompt photon production), by extending 
the model to other hadrons like kaons and hyperons, and relating their structure 
functions to fragmentation functions. This could also settle the issue about the 
parameters a and b. 
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Figure 1: (a) Comparison of the valence structure function data from ref.|p with the 
valence distribution within the statistical model (solid line). Dashed and dash-dotted 
lines are the valence densities obtained by the SMRS (ref.p) and GRS (ref.p) pa- 
rameterizations respectively, (b) The first two moments of the pion valence distribu- 
tion as predicted by the statistical model (solid line). Hashed region is the prediction 
[|] from the fit to Drell-Yan data of NAlO and the filled circles at = 49 GeV^ 



are lattice QCD |11] predictions. 
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Figure 2: Drell-Yan data for 7r~W ^ ix^ ii~X at beam momenta of 194 (left column) 
and 286 GeV/c (right column) taken by the NAlO collaboration 0. Solid lines are 
the next-to-leading order cross section predictions within the statistical model with 
the finite-size correction parameters a = —0.376 and b = 0.504 as used throughout 
this work. The dashed (a = —0.2, b = 0.3) and dot-dashed lines (a = —0.6, b = 0.7) 
at 286 GeV/c show the variation of the calculated cross sections with different choices 
of the parameters a and b. The K' factors for the 194 and 286 GeV/c data are 1.02 
and 1.06 respectively. 
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Figure 3: Drell-Yan cross section data for 7t~W —>■ fi^fi^X at beam momentum 252 
GeV/c taken by the E615 collaboration Solid lines are the next-to-leading order 
cross section predictions within the statistical model with a K' factor of 0.85. 
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Figure 4: Comparison of the statistical model parton distributions (solid lines) with 
those due to the SMRS (ref.p, dashed lines) and GRS (ref. 0, dash-dotted lines) 
parameterizations at = 30 GeV^. 
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